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Abstract. 

SciBooNE is a neutrino and anti-neutrino cross-section experiment at Fermilab, USA. The SciBooNE experiment is 
summarised and two independent CCQE analyses are described. For one of the analyses, an absolute -CCQE cross section 
in the neutrino energy region (0.6-1.6) GeV is shown and the technique developed for such a purpose is also explained. The 
total cross section measured over this energy range agrees well with expectations, based on the NEUT event generator and 
using a value of 1.21 GeV for the CCQE axial mass. 
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THE SCIBOONE EXPERIMENT 

SciBooNE is a and scattering experiment based at the Fermi National Accelerator Laboratory, USA. The main 
goal is to accurately measure the and cross-sections on nuclear targets at energies around 1 GeV. Currently 
few measurements have been made in this region on a nuclear target with high statistics. These measurements are of 
interest to the next generation of oscillation experiments. 



The Booster Neutrino Beam-line 

Protons are accelerated up to 8 GeV and focused onto a 71 cm long, 1 cm wide Beryllium target. The resultant 
meson shower (mostly pions) is then focused by a magnetic horn. The horn has reversible polarity allowing us to sign 
select the mesons that enter the 50 m decay pipe downstream from the target, this allows us to run in neutrino or anti- 
neutrino mode. The resultant mesons decay in flight producing our neutrino beam. The short decay pipe minimizes the 
chance of secondary decays resulting in an intrinsic v e background, important to oscillation studies. 



The SciBooNE Detector 

The SciBooNE detector is composed of three sub-detectors: the SciBar detector, acting as a target, the electromag- 
netic calorimeter (EC) and the muon range detector (MRD). 

The SciBar detector consists of 14,336 1.3 x 2.5 x 300.0 cm scintillator channels arranged in 64 layers comprised 
of one horizontal and one vertical plane read out by multi-Anode PMTs (MAPMT), with 64 channels read out per 
MAPMT. The SciBar detector was originally used by the K2K collaboration as a part of a near detector before being 
dis-assembled and shipped to Fermilab where it was reconstructed for use by SciBooNE. 

The EC is the central detector sandwiched between SciBar and the MRD. It consists of 2 planes - one horizontal, 
one vertical - each made up of 32 lead foil and scintillating fibre "spaghetti" calorimeter modules. Because of this, the 
EC is 1 1 Xq thick, able to measure large EM shower deposits exiting SciBar and thus identify n° and v e candidate 
events. The EC was originally built for the CHORUS experiment before being used by K2K. 



The MRD is a new detector made up largely of recycled materials. The detector consists of 362 scintillator paddles 
arranged into 13 alternating horizontal and vertical planes able to stop muons with momentum up to 1.2 GeV. The role 
of the MRD is to range out and reconstruct muons resulting from Charged-Current interactions. 

For the recycling effort and reuse of multiple detectors and electronic components the SciBooNE experiment was 
awarded a 2008 DOE P2* Environmental award. 



SciBooNE Operations 

SciBooNE ran from June 2007 until August 2008 collecting 0.99 x 10 20 Protons on Target (POT) in neutrino mode 
and 1 .53 x 10 20 POT in anti-neutrino mode, 25% more than requested. The total POT and event rate for the entire run 
can be seen in figureQ] 




FIGURE 1. The protons on target (POT) received during the full SciBooNE run (left) and the neutrino event rate per 4 x 10 
POT (right). The event rates during the two anti-neutrino runs are consistent. 



SciBooNE MC Simulation 

The neutrino flux MC is estimated using extrapolated thin target data from the HARP experiment, with the 
errors estimated using the Sanford-Wang parameterisation [1]. The NEUT library I2J is used to simulate the neutrino 
interactions with the nucleus in our MC simulation. Such a library is based on the following neutrino interaction 
models: The Llewellyn-Smith Model is used to model the CCQE interactions with an Ma value of 1.2 GeV/c 2 . The 
Rein-Sehgal model is used to model the CCln interactions, with the CC-coherent interactions suppressed by a factor 
of a third, in agreement with the SciBooNE measurement ft^. We use the Smith-Moniz Relativistic Fermi Gas nuclear 
model and Formation Zone Parameterisation for the final state interactions. For carbon the Fermi momentum is set to 
217 MeV and a binding energy of 27 MeV is applied. 



THE CCQE ANALYSES 
SciBar contained analysis 

Pre-selection cuts are used to identify SciBar contained events. The events must be completely contained within 
the SciBar fiducial volume (—130 < x,y < 130 and 2.62 < z < 157.2 cm) and the earliest track in the event must be 
within the beam window of 2fis. Charged-Current events are identified by tagging the decay electron signature from 
the stopped muon. This is identified by looking for double-coincidences between horizontal and vertical SciBar TDC 
blocks. Figure [2] shows the reconstructed muon lifetime in data for all identified events in SciBar. 

To purify the CCQE event fraction further cuts are applied. For the 1 -track sample a vertex activity cut of 8 MeV 
is applied to enhance the CCQE fraction with unreconstructed proton tracks, see figure [3] For multi-track events a 
Muon Confidence Level (MuCL) is constructed. Such a variable defines the probability of a particle being a MIP (or 
a muon-like track). A cut of 0.05 is applied to the track with the lowest MuCL value (least muon-like track), identical 
to the SciBar-MRD analysis (see figure |5). Those events that fail this cut become the 2-track fj, + % sample used in 
the fitting architecture to constrain the background. An additional CCQE cut is applied to the 2-track jj.+ p sample 
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FIGURE 2. Reconstructed muon lifetime in data. Lifetime: 2.003 ±0.047 (stat)[is consistent with experiment |4] 



TABLE 1. Summary table for the three SciBar 
contained sub-samples. 



Sample 


Data 


CCQE 


CCtt 


Other 


l" 


3733 


55.5% 


33.2% 


11.3% 


H +P 


349 


86.3% 


12.6% 


1.1% 


H + x 


1145 


34.0% 


47.9% 


18.1% 



(see figure |3]l using the reconstructed neutrino energy assuming a CCQE interaction and the energy deposited in the 
detector. The CCQE cut is defined as: 

{e ccqe + E ^ ScWar y Tp = £ (1) 

Where Ey C @ E is the reconstructed neutrino energy assuming a CCQE interaction, T p and are the proton kinetic 
energy and the total muon energy reconstructed from the energy deposits in SciBar respectively, and e is the unit-less 
parameter used for the cut. For this analysis a cut of e > —0.6 is applied to the sample. 




FIGURE 3. From left to right: The Vertex Activity cut applied to the 1-track sample, a cut of 8 MeV is applied. The CCQE cut 
applied to the 2-track fl + p sample, a cut of -0.6 is applied. 

The final sample selection is summarised in table Q] The 1-track sample is the largest statistical sample, with the 
2-track (/x + p) having a high CCQE purity and the 2-track (jj. + n) a high background purity. Figure H shows the 
Pjj and cos(0|,) distributions for all three samples. The 2-track sample contains backwards tracks in data and MC, 
something under investigation. 
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FIGURE 4. Pp (top) and 0^ (bottom) distributions (CCQE - white, background - hashed, dirt events - black) for all 3 final SciBar- 
contained samples, from left to right: 1-track [A events, 2-track /x + p events, 2-track /x + it events. Dirt events occur when particles 
from a neutrino interaction outside the detector are detected within the detector. This is a minor background in the SciBar-contained 
analysis and is negligible for the SciBar-MRD analysis. 



SciBar-MRD analysis 

Analogous to the previous analysis, the events in this analysis must satisfy two pre-selection cuts. The events must 
be produced within the beam time window (0 < t < 2/is) and the interaction vertexe^] must be inside of the SciBar 
fiducial volume (—130 < x,y < 130 and 2.62 < z < 157.2 cm). The main difference to the SciBar-contained analysis 
is the requirement that the events contain at least one SciBar-track reaching the MRD detector (SciBar-MRD matched 
sample). Those tracks are defined as muon-track candidates. In addition, the muon tracks are required to stop inside of 
the MRD detector in order to reconstruct completely the neutrino energy, obtained under the assumption of a CCQE 
interaction. The resultant sample is called the SciBar-MRD stopped sample. 

The CCQE interactions are characterised by the low track multiplicity, the track vertex connection and the presence 
of a proton track. These three properties determine the CCQE cut selection. The track multiplicity cut (fig. [5} takes 
advantage of the QE event topology, selecting events with two tracks. However, the 1-track events are selected as well 
in order to include the QE events with low energy proton-tracks not reconstructed (which correspond with a proton 
momentum threshold of around 450 MeV/c). 
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FIGURE 5. From left to right, the distribution of number of tracks for the SciBar-MRD stop sample, the distribution of the 
connection vertex for the 2-track sample and the distribution of the MuCL variable for the 2-track sample after passing the 
connection vertex cut. The Monte Carlo is normalised to the SciBar-MRD matched sample. 

The connection vertex and the particle identification cuts are based on properties of the second track, then both are 
applied to the 2-track sample. The connection vertex is defined as the distance between the neutrino interaction vertex, 
defined above, and the upstream edge of the second track. For QE interactions, this distance must be close to zero. The 
figure [5] shows the connection vertex distribution for the 2-track sample with cut value of 10 cm. This value has been 
selected by maximizing the efficiency times purity of the cut. 



The interaction vertex is defined as the upstream edge of the muon-track candidate 



TABLE 2. Summary table with the three sub-samples selected for the 
extraction of the CCQE cross section. 



data CCQE CCQE ratio 

Sample events Efficiency (%£ Purity (%) data/mc[_ 

1 track 13586 52.9 65.2 1.03 

2track(ji+p) 2915 11.1 68.5 1.05 

2 track (n + n) 1628 2.7 32.3 1.14 



* calculated with respect to the number of CCQE events in the SciBar-MRD 
matched sample 

1 The monte carlo is normalized to the SciBar-MRD matched sample. 



The particle identification cut classifies the second track as proton-like or pion-like track according to a deposition 
energy study. For such a purpose, the MuCL variable (defined in the previous section) is used. Figure [5] shows the 
MuCL distribution of the second track coming from 2-track events which pass the connection vertex cut. This figure 
shows that second tracks with MuCL values less than 0.05 are defined as proton-candidates while higher values 
correspond to pion-candidates. 

The final selection is composed of three samples (see table [2j. The 1 -track and 2-track (/I + p) samples define the 
QE part of the selection and are characterised by a high statistics and high CCQE purity respectively. The 2-track (ji 
+ n) sample defines the nonQE-part of the selection. This sample, mainly composed of CC resonant K + production, is 
used in the fit to constrain the overall amount of CC nonQE background interactions with data, as explained below in 
the fit method description section. 



EXTRACTION OF ABSOLUTE CCQE CROSS SECTION 
Fit method description 

The fit method is based on minimizing a likelihood function with the assumption that observed events follow a 
Poisson distribution. The likelihood expression is then written as follows |@]: 

X 2 = 2 x EEC -N» +N/f s x In , (2) 

where Nff p and N°j correspond to the expected and the observed number of events represented in bins of reconstructed 

muon momentum (p") and reconstructed muon angle (9^). The number of expected MC events is mathematically 
expressed as follows: 

k 

that is, the MC is broken into a QE component (Nfj E ' k ) which also depends of the true energy bins (given by the index 

k) and a non-QE component (N"""® E ). By taking expression ©, the fit parameters are introduced in the following 
way: 



*% P = F N 



i=10 

L\jQE.k . KT nonQE 
a k Nfj +a bck N ij 

k 



(4) 



A total of twelve parameters define the fit, two of them are global parameters (Fn and a^) and the other ten are 

neutrino energy dependent parameters (fli ,«2 ,«io)- Fn is a scale factor, containing information about the data/MC 

normalization, gives information about the non-QE contamination. Finally, the energy dependent parameters (a k ), 
each one associated with a different true neutrino energy bin, return information of the neutrino energy shape. 

The absolute CCQE cross section is then evaluated for each neutrino energy bin using the following formula: 



Ov„ [E k v )=F N xa k x a^ EUT (E v ) , (5) 

where Oy EUT (Ey) corresponds to the value of the absolute cross section predicted by NEUT in the neutrino energy 
bin k. One can observe again the role of the a% parameters in the expression (0, re-weighting the predicted cross 
section for each neutrino energy bin. 



Results 

The absolute v^-CCQE cross section per nucleon shown in this proceeding (see fig. [6j is extracted only by using 
the sample selection from the SciBar-MRD analysis. The neutrino energy range of the cross section is then limited by 
the sensitivity of this selection, which only works in the interval (0.6-1.6) GeV. The low energy region will be covered 
by adding the SciBar-contained sample (work in progress). 

The cross section results are consistent with the theoretical prediction given by NEUT[2], as figure [6] shows. The 
systematic uncertainty is due to the limited knowledge in the neutrino flux predictions at SciBooNE. The dominant 
contribution to the flux systematics comes from the uncertainty in the rate and kinematics of 7C + -production in 
the interactions of beam protons with the beryllium target, as measured by HARP [6],|1]. Other flux systematics 
considered include uncertainties in other hadronic interactions, in beam proton delivery, and in horn magnetic focusing 
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FIGURE 6. The absolute v^-CCQE cross section measured in SciBooNE for the SciBar-MRD analysis (circle dots) in compar- 
ison with the NEUT prediction (triangle dots). The statistics (bars) and the flux systematic errors (grey band) are included in the 
data. 
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